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In a three-dimensional simulation higher-order derivative correlations, 
including skewness and flatness (or kurtosis) factors, are calculated for 
velocity and passive scalar fields and are compared with structures in the 
flow. Up to 128 3 grid points are used with periodic boundary conditions 


in all three directions to achieve R\ to 82.9. The equations are forced 
to maintain steady-state turbulence and collect statistics. The scalar- 
derivative flatness is found to increase much faster with Reynolds number 
than the velocity derivative flatness, and the velocity- and mixed-derivative 
skewnesses do not increase with Reynolds number. Separate exponents 
are found for the various fourth-order velocity-derivative correlations, with 
the vorticity-flatness exponent the largest. This does not support a major 
assumption of the lognormal and models, but is consistent with some 
aspects of structural models of the small scales. Three-dimensional graphics 
show strong alignment between the vorticity, rate-of-strain, and scalar- 
gradient fields. The vorticity is concentrated in tubes with the scalar 
gradient and the largest principal rate of strain aligned perpendicular to 
the tubes. Velocity spectra, in Kolmogorov variables, collapse to a single 
curve and a short -5/3 spectral regime is observed. 

§A. Introduction The classical approach to investigating small-scale intermit- 
tency in turbulence is through the higher-order derivative correlations such 
as skewness and flatness factors. Experimentally this has been done for both 
the velocity and temperature, which is a passive scalar when buoyancy is 
negligible. The long-range goal is to improve our understanding of the struc- 
ture of the small scales, possibly leading to improved methods for subgrid 
modelling. An intermediate objective has been to relate the derivative cor- 
relations to dissipation correlations and corrections to Kolmogorov scaling. It 
is believed that this is possible because the small scales are universal; that is, 
the small scales have a structure that is independent of the large scales and 
can be modeled. So far only single-probe velocity and scalar statistics, such as 
the velocity-derivative flatness and skewness, have been measured with hot- 
wire instruments. More complicated statistics have not been measured and 
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alternative approaches, such as flow visualization, are limited because they 
do not have the flexibility necessary to distinguish small-scale structures that 
are intermittent in space and time. 

Another approach to investigating the small scales is numerical simulation. 
In a simulation more detail can be obtained than in experiments and the con- 
ditions can be more closely controlled. For example, because all components 
of the velocity are known, one is able to study correlations beyond the deriva- 
tive skewness and flatness and computer graphics can display structures not 
accessible to experiments. The main disadvantage of a simulation is that only 
a limited range of length scales is allowed, which restricts the Reynolds num- 
ber to very low values. But in low-Reynolds-number experiments the values 
of most of the derivative correlations are significantly different from their 
uncorrelated or Gaussian values. Therefore, these statistics are accessible to 
current numerical methods and computers. Siggia (1981a) used a numerical 
simulation to calculate higher-order correlations of the velocity, and discussed 
the relation between small-scale vortex structures and intermittency. Our 
approach will be to use a similiar simulation to look in more detail at the 
velocity-derivative statistics and the statistics of a passive scalar. Graphical 
display of small-scale vorticity, rate-of-strain, and scalar-gradient structures 
is used to interpret these statistics and comparisons with phenomenological 
theories and experiments are made. 

There are two phenomenological approaches to predicting the small-scale 
statistics: either by assuming a form for vortical structures or a form for the 
energy cascade from large to small scales. Two models that are based on the 
cascade of energy are the lognormal theory of Kolmogorov (1962) and the 
/?-model of Frisch, Sulem, and Nelkin (1979). Both theories predict a cor- 
rection to the A; -5 / 3 inertial-range kinetic-energy spectrum of Kolmogorov 
(1941). They also predict that as the Reynolds number grows that the 
velocity fluctuations become increasingly localized, or intermittent, distribu- 
tions become highly non-Gaussian, and the higher-order correlations, such 
as the derivative skewness and flatness factors, increase with Reynolds num- 
ber with a power- law dependence. The power-law exponents depend on the 
details of each model and on fi, the characteristic exponent of the dissipation- 
dissipation correlation function in the inertial subrange, rj < < r < < L, 

< e(x)e{x + r) >=Ae 2 {L/rY (1) 

(Monin and Yaglom 1975, p. 618), where A is a constant and L is defined 
by (12). All correlations of a given order are predicted to have the same 
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power-law dependence. The lognormal model predicts that 


< € m (x)e n (x + r) >~ (L/r) Mmn 


( 2 ) 


and that 



where 

3 

= -t*n{n — 1) (4) 

4 

(Frenkiel and Klebanoff 1975). The /?-model is similiar, but predicts that 


Ct n 


3 ju ( n 
4-/T2 


1) 


( 5 ) 


(Nelkin and Bell, 1978). The success of these models in predicting the ex- 
perimentally observed dependence of derivative skewness and flatness on 
Reynolds number is discussed in detail by Antonia, Satyaprakash, and 
Hussain (1982), who conclude that the lognormal model is superior in this 
respect. Neither of the corrections to the A : -5 / 3 law proposed has been ob- 
served, although the lognormal correction is much smaller and probably not 
much above statistical noise. 

The value of fi can be found either by direct measurements or by using these 
models to calculate back from a n . Nelkin (1981) and Antonia, Satyaprakash, 
and Hussain (1982) summarize the current evidence and conclude that most 
previous estimates were too high. They suggest that n = 0.2. 

Structural models assume that the small scales are composed of tubes 
or sheets of vorticity. Corrsin (1962) assumed that sheets dominated and 
concluded that <24 = 1.5. Saffman (1968) also assumed that the vorticity 
would be found in sheets whose thickness would be the order of the Taylor 
microscale (9), but that within the sheets, dissipation would be localized in 
regions characterized by the Kolmogorov length scale (10a). He found that 
o 4 = 1 and o 3 = 0, that is, the skewness is constant. A variation, proposed 
by Tennekes (1968), assumes that the dominant structures are tubes whose 
thickness is the Kolmogorov length scale, but which are subject only to the 
large-scale strain. His results agree with those of Saffman. Experimentally, 
o 3 is observed to be very small, possibly zero, but a 4 is much less than the 
predicted value of one. 
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Dissipation of the variance of a passive scalar is observed to be intermit- 
tent, much like kinetic-energy dissipation. For example, the scalar-derivative 
flatness scales with Reynolds number in a manner similiar to the velocity- 
derivative flatness, but with a larger exponent (Antonia and Chambers 1980). 
Sudden jumps in experimental temperature signals, known as “ramps”, are 
also observed (Antonia et al. 1979). The only theoretical attempt to describe 
scalar intermittency was by Van Atta (1974), who used lognormal assump- 
tions. Correlations between the velocity and scalar derivatives would also be 
interesting, but because it is difficult to measure two turbulent fields simul- 
taneously, there have been almost no experimental measurements of these 
correlations. 

In this paper, a variety of velocity- and scalar-derivative correlations are 
calculated over a range of Reynolds numbers. In an isotropic uncorrelated 
field with Gaussian statistics each of the correlations to be discussed has an 
easily determined value. For the velocity-derivative correlations Fq u » (3), 
the skewness (n=3) is zero, the flatness (n=4) is 3, and the sixth order 
correlation is 15. For mixed correlations in which both the velocity derivative 
and the scalar derivative have a power of 2, such as the mixed-derivative 
correlation (40), the uncorrelated value is 1. At high Reynolds numbers the 
statistics are usually highly non-Gaussian, both in experiments and in our 
simulations. We will concentrate on third- and fourth-order correlations, with 
some fifth- and sixth-order correlations presented to allow comparisons with 
the phenomenological models. Since the models discussed predict scaling laws 
and exponents, we will estimate these exponents and make comparisons. 

The governing equations are the incompressible Navier-Stokes equation for 
the velocity and the transport equation for a passive scalar. The Navier- 
Stokes equation is 


The nonlinear term, u • Vu, can be written in several different forms which 
are computationally convenient: the conservative form, 


du/dt -f - u-Vu = — Vp/p i'V 2 u 


(6a) 


V • u = 0 (incompressibility) 


V • ( uu ) 


(6b) 


and the rotational form 


u X u + V£u 2 


(6c) 
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The scalar equation can also be written two ways in incompressible flow: the 
convective form 

dti/dt + u- V0 = DV 2 ® (7a) 

and the conservative form 

dti/dt + V • u& = DV 2 & (7b) 

In the absence of viscosity v and diffusivity D the equations conserve 
two positive-definite quadratic invariants: the kinetic energy of turbulent 
fluctuations, 

E = - < tijUj > 

and the scalar variance, 

E e =< 6 2 > 

where the domain of integration is over all space. The fundamental dimen- 
sionless parameters that we use are the Taylor-microscale Reynolds number, 
R\ = UX/z/, and the Prandtl number u/D, where U is the characteristic 
velocity of the turbulence, 


-U 2 = E 
2 


( 8 ) 


and X is the Taylor microscale, 


( 9 ) 


X = < «i 2 > 1 ^ 2 / < (dui/dxi ) 2 > x t 2 

Also of interest are the kinetic-energy dissipation rate, 

Cl Q rrO 

e = ~It < $ U > 

and the scalar-variance dissipation rate, 

d .o 

x = ~jt <e > 

The Kolmogorov length scale, wavenumber cutoff, and velocity scale are 
defined as 

»?=(—)*, K k = ~, and ( a ') 1 / 2 (10a) 

e V 

and the Oboukov-Corrsin microscale and wavenumber (Corrsin 1951) are 
defined as 

Voc = (—)i and Koo=— (10b) 

€ VOC 
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§B. Numerical Method 

Two numerical codes were used for the simulations presented. Each is 
a three-dimensional spectral code with periodic boundary conditions. By 
spectral we mean that the fundamental variables that are stored and ad- 
vanced in time are the Fourier-transformed velocity and scalar fields, u(k ) 
and 9(k). For simple geometries these methods allow more resolution than 
finite-difference methods for the same number of grid points (Orszag 1971). 
To improve the speed of spectral methods fast Fourier transforms are used to 
return to physical space, where the nonlinear terms are calculated by forming 
products. This method is sometimes referred to as pseudospectral because it 
introduces aliasing errors in the high wavenumbers. One way to reduce these 
errors, or dealias, is by truncating interactions outside a boundary in Fourier 
space. This does not significantly affect the resolution because only the high- 
wavenumber tail of the spectrum, with only a small fraction of the total 
energy, is truncated. To completely dealias one should truncate all wavenum- 
bers larger than N / 3, where N is the number of grid points in a single direc- 
tion. To completely dealias a three-dimensional code without truncation re- 
quires calculating the nonlinear terms on eight shifted grids. Fortunately, 
Patterson and Orszag (1971) have shown that in two or more dimensions all 
the aliasing errors can be eliminated by using only two shifted grids and a 
spherical truncation with the proper radius. This method was used in their 
original code, known as Superbox. Another approach uses shifted grids on 
alternate evaluations (Rogallo, 1981). The extent to which aliasing errors 
affect a calculation also depends on the algorithm used for the convective 
terms. Superbox uses the aliased rotation algorithm (6c) because it conserves 
kinetic energy, which helps inhibit instabilities and reduces the need for using 
shifted grids. Comparisons between versions of Superbox with and without 
shifted grids show no significant differences in their spectra (Patterson 1980). 

The scalar equations do not allow aliasing to be neglected as easily as 
the Navier-Stokes equations do. For zero diffusivity, the scalar variance 
is conserved by the exact equations. But the aliased versions of both the 
convective and conservative equations (7a, b) do not conserve scalar variance. 
Instead, their respective aliasing errors are equal in magnitude and opposite 
in sign. We investigated taking advantage of this by averaging the two 
forms of the scalar equation to dealias, but it did not yield any significant 
improvements because our fields were already well resolved. Arakawa (1962) 
used similiar considerations in writing his stable finite-difference scheme for 
the two-dimensional Navier-Stokes equations. 
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The code used for the 32 3 runs is based on the Superbox code of Siggia 
and Patterson (1978) and Siggia (1981a). The rotational form of the Navier- 
Stokes equation (6c) was used for the velocity and the time advancement was 
leapfrog with stabilization done every 40 time-steps by a second-order Runge- 
Kutta step. Shifted grids were not used for dealiasing, but the convective and 
conservative forms of the scalar equation (7a, b) were used on alternate time- 
steps to partially dealias the scalar. Wavenumbers were truncated outside a 
sphere of radius N/ 2. This truncation has the advantage of being isotropic, 
but does allow significant one-dimensional aliasing errors. Details of this 
code may be found in Kerr (1981). Simulation of the velocity and one scalar 
required approximately 0.75 second of computer time (cpu) per evaluation on 
a Cray-1 S computer and the longest simulation used 40 minutes of computer 
time. 

The code for the 64 s and 128 s calculations also used the rotational form of 
the Navier-Stokes equations, with a spherical truncation of N / 2 and without 
shifted grids, but only the conservative form of the scalar equation was used. 
The time advancement was third-order Runge-Kutta (Wray 1981) and up to 
three scalars could be calculated simultaneously. Two-to-one word packing 
was used to improve data management (or I/O). For a 64 3 mesh with three 
scalars, 6.5 seconds of Cray cpu time and 2.0 seconds of I/O wait (time waiting 
for retrieval of data from disks) were required each evaluation, with three 
evaluations per time-step. The code with no scalars required 3.25 seconds of 
cpu time and 1 second of I/O wait per evaluation. Our longest simulation for 
a 64 3 mesh with three scalars (F19-21) required 12 hours of computer time 
to simulate seven eddy-turnover times (11). 

We tried to run a 128 3 mesh using Fortran I/O on the Cray-lS with disks, 
but due to the large number of files on disk, I/O wait became intolerably large. 
What made possible a 128 3 simulation was the installation of a solid-state 
storage device (SSD), which effectively expanded the memory of the Cray-lS 
to 17 million words. For a 128 3 mesh with three scalars approximately 55 
seconds of cpu time and a negligible amount of time waiting for I/O between 
the SSD and main memory were required per evaluation. Three eddy-turnover 
times required 70 hours of Cray time. 

§C. Forcing, Simulation Parameters, and Aliasing Limits 

The codes outlined have been used for forced and decaying simulations. 
This paper will concentrate on the higher-order correlations obtained when 
the large scales are forced. The decaying calculations are discussed in detail 
by Herring and Kerr (1982) and Kerr (1981), though some of those results 
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are included here in figure 3 and table 3. 

Why do we force the large scales? Ideally we would like to simulate an 
infinite domain with unlimited resolution of the smallest scales, but in prac- 
tice we are limited to a finite box with a finite mesh. There is a minimum 
wavenumber associated with the size of the box and a maximum wavenum- 
ber associated with the mesh spacing. If the effects of these limits are to 
minimized, the energy must be restricted to a very small range of wavenum- 
bers and the Reynolds number must be small. The calculations discussed 
in Herring and Kerr (1982) simulate decaying flow behind a grid and are 
restricted in this way. 

By forcing the large scales, some of the restrictions imposed by the low- 
wavenumber limit can be removed. It is believed that for sufficiently large 
Reynolds number that there is a cascade of energy from large to small scales 
which maintains a statistically-steady state. By forcing we hope to mimic a 
cascade from wavenumbers smaller than the minimum wavenumber to those 
in our computational box. Forcing would also maintain a steady spectrum at 
moderate Reynolds numbers and allow us to collect statistics of the higher- 
order correlations. Since the derivative correlations are representative of the 
small scales or high wavenumbers, it is hoped that they will not be affected by 
the details of the large scales, in our case the details of the artificial forcing 
used. Our approach is opposite to that of many turbulence modelers, who 
are primarily interested in the large scales and model the small scales. 

Both the velocity and scalar were forced by time advancing the lowest 
wavenumber band ( 1 < k/k 0 < 2, where k 0 is the lowest wavenumber of 
the code) independently of the high wavenumbers. The modes calculated in 
this manner were included in the full equations and the higher modes were 
advanced as usual. Details of the numerics may be found in Kerr (1981). 
Analysis of single modes in the forced wavenumber band indicates that the 
modes behave chaotically, although a complicated periodic motion cannot be 
ruled out. Our forcing is different than that used by Siggia (1981a), whose 
forcing modes were outside his computational box. 

We present results from 27 forced cases for the Taylor-microscale Reynolds 
number, R\, between 9 and 82.9 and the Prandtl number between 0.1 and 
2.0 (see table 1), but will discuss only a few in detail. All but one of the 
simulations includes a passive scalar. The exception (labeled F12 in table 1) 
uses the original dealiased Superbox code and our forcing to discern the effect 
of aliasing. 

Each simulation was run until a statistically-steady state was reached, 
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usually at 
as 


about two eddy-turnover times. An eddy-turnover time is defined 


te = 


L 

U 


( 11 ) 


where 

L= [ k~ l E(k)dk 
4E J 


( 12 ) 


is the integral length scale. After reaching a steady state, statistics were 
collected for at least one and a half more eddy-turnover times and the means 
and variances of a variety of velocity and scalar correlations were obtained by 
averaging. These statistics and the errors are discussed in detail in appendix 
A. 

Most earlier 32 s spectral simulations (Orszag and Patterson 1972) have 
been limited to R\ less than 30 based on considerations of the location of the 
peak of the three-dimensional velocity-dissipation spectrum (which is near the 
“bump” discussed in the following section for figure 1). Further restrictions 
come from considerations of trends in the derivative flatnesses. Generally 
the derivative skewnesses (14,15) and flatnesses (32,39) increase, or remain 
constant, with increasing Reynolds number. Failure to do so indicates that 
small-scale truncation and aliasing errors are significant. This behavior was 
observed for a 32 3 mesh at Reynolds numbers larger than 24. In order to 
demonstrate the magnitude of aliasing errors at this Reynolds number, table 
2A gives derivative flatness factors for 32 s and 64 3 calculations starting with 
identical initial conditions at ^29. On the 64 3 grid there are virtually no 
aliasing or truncation errors at this Reynolds number for any of the Prandtl 
numbers used. Any differences between the 32 s and 64 s results are due to 
aliasing and high-wavenumber truncation errors in the 32 3 simulation. In 
run Fll only spherical truncation is used to dealias; while F12 uses the 
original Superbox code to dealias by shifting grids on each evaluation. For 
Fll(aliased) both the velocity-derivative flatness and scalar-derivative flatness 
for Pr = 0.5 are significantly lower than the 64 3 results (F14), while in 
F12(dealiased) there is virtually no difference with FI 4. Therefore, at this 
Reynolds number an aliased 32 s calculation does not resolve the small scales, 
whereas a dealiased calculation does. 

To estimate the maximum allowable Reynolds numbers on the 64 3 and 
128 3 grids, the Kolmogorov wavenumber cutoffs of the simulations can be 
compared to the maximum wavenumbers allowed by the grids. For the inertial 
range ( E(k ) ~ K 0 € 2 / 3 k~ 5 / 3 ) R\ ~ K ^ . But the Reynolds numbers in 
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the 32 3 and 64 3 simulations are too low for an inertial range to appear and 
R\ ~ Kk 1-25 fits our results better. If the maximum allowable Kolmogorov 
wavenumber on a 64 3 grid is assumed to be twice that permitted on a 32 3 
grid, then to ensure that the small scales in the 64 3 simulation are fully 
resolved the maximum R\ is 55.9(~ 24 X 2 1 - 25 ). As shown in figure 1, our 
128 3 simulation has a noticeable -5/3 inertial range, so the K ** law should 
be used and f? x = 82.9(~ 52 x 2°- 67 ). 

In addition to its effect on the calculation of the velocity and scalar fields, 
aliasing can introduce errors in the calculation of the derivative flatness 
factors. To estimate this error the derivative flatnesses for a 64 3 simulation 
were calculated on a 32 3 grid for one realization. Table 2B shows that the 
velocity-derivative flatness calculated on a 32 3 grid yields a slightly smaller 
result. This reinforces our conclusion that R\ 29 is slightly too high 
for a 32 3 grid. The scalar-derivative flatness for Pr — 0.5 shows aliasing 
errors similiar to those in the velocity field, so no further restrictions on 
the Reynolds number are necessary when Pr < 0.5. On the other hand, a 
significant decrease is seen in the scalar-derivative flatness for Pr = 1.0 when 
it is calculated on a 32 s grid. 

Naively one would expect that for Prandtl number of 1, where the vis- 
cosity and scalar diffusivity are the same, that there would be similiar trun- 
cation effects for the velocity and the scalar. Instead, the scalar errors are 
much larger. What does theory tell us? For Prandtl number greater than 
1, Batchelor (1959) predicts a new spectral regime in which scalar variance 
diminishes as k~ 1 for wavenumbers greater than the Kolmogorov wavenum- 
ber. The cutoff for the new spectral regime is ks = (e/uD 2 ) 1 ^ 4 . This 
spectral regime has been observed (see Monin and Yaglom 1975, p.513). Since 
a higher Prandtl number implies a higher spectral cutoff, for large Prandtl 
numbers truncation effects will appear even at very low Reynolds numbers. 
Consequently for Pr = 2, R\ > 12 was not run on a 32 3 mesh. But an 
anomalous k~ l subrange is also observed experimentally for Prandtl number 
as low as 0.7, the atmospheric value (the “bump”, Hill 1978). Therefore, even 
at Pr = 1 truncation, as indicated by the scalar-derivative flatness, becomes 
important for R\ > 18 on a 32 s grid. 

Because of these difficulties our simulations are restricted to low Prandtl 
numbers. Some of the simulations for Pr = 1 are included in the plots of 
the correlations, but only for lower Reynolds numbers. Prandtl numbers of 
current experimental interest that we can simulate are 0.7, the atmospheric 
value, and moderately low Prandtl numbers, such as found in liquid metals. 
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Extremely low values cannot be simulated satisfactorily because the relevant 
low-wavenumber cutoff is the Oboukov-Corrsin cutoff, Koc (10b), and this 
would be below the minimum wavenumber of the simulation.The lowest value 
we use is Pr = 0.1. Analysis of spectra in the next section support these 
conclusions. 

§D. Spectra and Skewnesses 

Since most simulators determine whether there is significant truncation 
error by inspecting high-wavenumber spectra, we have included some spectra 
to demonstrate the validity of our method. Figure 1 shows three-dimensional 
energy spectra, normalized by the Kolmogorov microscales and multiplied 
by A; 5 / 3 , for several Reynolds numbers. The turnup at large wavenumbers 
is a numerical effect of the large-wavenumber cutoff. Where the turnup is 
small, the small scales are well resolved. The turnup for R\ = 55.9 in a 
64 s simulation and R\ = 82.9 in the 1 28 s simulation is the most we could 
accept. The three-dimensional kinetic-energy and scalar-variance spectra for 
our largest Reynolds number, 82.9, are plotted in figure 2. The velocity 
spectrum is normalized as in figure 1 and the scalar spectra are normalized by 
the Kolmogorov microscales and the scalar-variance dissipation in a similiar 
manner. When compared with the scalar variance at low wavenumbers, the 
turnup at high wavenumbers for Pr = 1 is unacceptably large. This reinforces 
our conclusion that moderate-to-high Prandtl numbers can be simulated only 
with low-Reynolds-number simulations. 

With the exception of deviations in the lowest Reynolds-number spectrum 
(R\= 18.5), all of the spectra in figure 1 collapse to a single curve in the 
dissipation range. In addition, a short -5/3 regime appears in the 128 3 
spectrum (R\=S 2.9) at low wavenumbers. For comparison, the dashed line 
in figure 1 fits a spectrum of the form predicted by Pao (1965), 

E{k) = K 0 6 2 ^k-^ s exp(-1.5K 0 ue-^ 3 k^ 3 ) (13) 

to the computed spectra with K 0 , the Kolmogorov constant, equal to 2.45. 
Although experimental values for the Kolmogorov constant as high as 2.1 
have been reported (Gibson et al. 1970), the fitted value is much higher 
than the usual experimental value of 1.4 to 1.7 (Monin and Yaglom 1975, 
p.483-485). To understand the apparent differences between the calculated 
and experimental Kolmogorov constant the velocity-derivative skewness must 
be considered. From this we will conclude that fitting the calculated spectra 
with Pao’s form could lead to anomalously high values for the Kolmogorov 
constant. 
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For both the velocity and scalar three classes of the derivative skewnesses 
can be defined. These are the real-space derivative skewnesses, spectral- 
transfer skewnesses, and dissipation skewnesses. The real-space derivative 
skewnesses are measurable with hot-wire probes. The spectral and dissipation 
skewnesses are not strictly third-order moments, but in isotropic flow they 
are related to the real-space derivative skewnesses. The measurable real-space 
derivative skewnesses are the velocity-derivative skewness, 



the mixed-derivative skewness 



and the scalar-derivative skewness 


< ( ffi ) 3 > 
<(^) 2 > } 


(16) 


In isotropic turbulence the terms in the velocity-derivative skewness and the 
mixed-derivative skewness can be written more generally. That is, 



4 o 

- < tr e 3 > = 
3 


(Betchov 1956) and 



9u i, 96 )2 e p Q 

dx l [ dx x ) 


> 


Furthermore, 

< (dui/dxxf >= ^e/u 


(17a) 

(17b) 


(18) 

(19a) 
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and 


(19b) 


<( 


89 

dxi 



These rotationally-invariant forms are used in all our calculations of the real- 
space skewnesses. 

The spectral-transfer skewnesses and the dissipation skewnesses are related 
to the real-space skewnesses through integral equations for the dissipation of 
energy and scalar variance, 


-r<=( 
2i/ dt J 

f k 2 T u (k)dk -2u J 

^ k 4 E u (k ) dk 

(20) 

1 d _ 

2 Ddt X ~ . 

j k 2 Tg(k)dk - Dj 

[ k 4 E$(k) dk 

(21) 


where T u (k) and Tg{k) are the nonlinear energy and scalar-variance trans- 
fers, respectively. For our forced equations another term should be added, 
but this term is negligible except at low Reynolds numbers. The second- 
order moments of the transfers are related to the real-space skewnesses by 
f k 2 T u (k)dk =< uiieijUj > and f k 2 Tg(k)dk = — < O^ijOj >, identically. 
The integrals of the transfers are zero (/ T(k)dk = 0), but because there is a 
net cascade of kinetic energy and scalar variance to higher wavenumbers, their 
second-order moments will in general be positive and nonzero. Therefore, we 
can define two spectral-transfer skewnesses, the velocity spectral skewness 


c _ 2 / k 2 T u (k)dk 
and the mixed spectral skewness 


2 Pn 

35 (T3F) 


S$(k) — 


2 Sj^Tg{k)dk 
15 (t£f)*(§£») 


( 22 ) 


(23) 


where Pq is the rate of enstrophy production and 0, the enstrophy, equals 
^ < UiUi >. Because the spectral-transfer skewnesses are nonzero, in 
isotropic flow the velocity- and mixed-derivative skewnesses will also be non- 
zero and obey 

S u - Su(k) Sud == ^8(k) 
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These arguments do not apply to the real-space scalar-derivative skewness 
(16), which in an isotropic calculation such as ours is zero. However, it has 
recently been found that in shear flows with a temperature gradient that the 
scalar-derivative skewness is nonzero (Sreenivasan and Tavoularis 1980). 

For our statistically-steady forced simulations -+■ 0 and ->• 0, so 
we expect that as the Reynolds numbers increases that f k 2 T u -*• v f k 4 E u 
and / k 2 Te -*• D f k 4 E$. Therefore, a third measure of the skewnesses, the 
dissipation skewnesses, can be defined by (Wyngaard 1971) 


S e = 


2 e n 

35 <Tfe)* 


(24) 


and 


where 


s _ 2 D [k^e(Jc)dk 
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at J 


(25) 

(26) 


is the enstrophy dissipation rate. 

As discussed in the introduction, all the higher-order correlations are ex- 
pected to diverge from their uncorrelated values as the Reynolds number in- 
creases. Experimentally, the velocity-derivative skewness increases from zero, 
the uncorrelated value, to a value greater than 0.3 at very low experimental 
Reynolds numbers and appears to increase indefinitely in very high Reynolds- 
numbers atmospheric measurements (Tavoularis et al. 1978). But at inter- 
mediate Reynolds numbers the experiments are less conclusive. Tavoularis, 
Bennett, and Corrsin (1978) cite data which suggests that the skewness 
decreases with Reynolds number between R\ — 10 and 500. Klebanoff (1982) 
finds the velocity-derivative skewness to be constant and equal to 0.43 between 
R\ = 40 and 200 in a boundary layer. This experimental result is plotted in 
figure 3. 

Simulations discussed in Herring and Kerr (1982) show that the velocity- 
derivative skewness increases from zero to 0.4 by R\ ~ 20, in agreement 
with low-Reynolds-number experiments. They find that the mixed-derivative 
skewness increases in a similiar manner. At intermediate Reynolds numbers, 
they suggest that both the velocity- and mixed-derivative skewnesses reach 
maximum values, then remain constant with increasing Reynolds number, 
agreeing with KlebanofFs experiment. With our 64 3 and 128 3 simulations we 
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have determined the skewnesses at higher Reynolds numbers and confirmed 
those results (figure 3). Figures 3 and 4 plot the real-space and dissipation 
skewnesses, respectively, as functions of R\ and in all cases the magnitudes 
of the skewnesses converge to approximately 0.5. The value for the velocity- 
derivative skewness for simulations Fll to F27 is —0.505 ± 0.005. 

The statistical models of the small scales and the atmospheric measure- 
ments mentioned predict that the velocity-derivative skewness will increase 
very slowly with Reynolds number (4,5). Although the accuracy of our curve 
is not good enough to rule out the existence of such a trend, it appears 
significant that both our simulation and some experiments show no Reynolds- 
number dependence over a wide range. The reason the velocity-derivative 
skewness does not increase with Reynolds number is probably related to the 
alignment of small-scale structures. Recall that the model of Tennekes (1968) 
suggested that if the structures are tubes, then the skewness would not in- 
crease with Reynolds number. In section G (Graphics) it will be seen that 
our structures are tubes. In addition, fourth-order statistics and our graphics 
find strong alignment between the vorticity and the rate of strain. 

Because of the connection between the skewnesses and the fourth-order 
moment of the spectra, a simple relation between the velocity-derivative 
skewness and the Kolmogorov constant exists. If one assumes Pao’s spectral 
form (13) and uses the equation for the velocity-dissipation skewness(24), then 

S u = —^ni^n(^K < ,)=? = -2A0K^ (27) 

For S u = —0.5, K 0 = 2.85. For a typical experimental Kolmogorov con- 
stant, K 0 = 1.7, (27) gives S u = —1.08, which disagrees with experimental 
measurements of the skewness at moderate Reynolds numbers. If the ex- 
perimental values for the Kolmogorov constant and the skewness are correct, 
then Pao’s spectral form is clearly a poor approximation to the experimental 
spectra. There are also significant differences between the calculated spectra 
of figure 1 and Pao’s form near ( rjk ) ~ 0.2. If the calculated Kolmogorov 
constant is to be consistent with the experiments (that is K 0 ~ 1.7), then 
at wavenumbers below those simulated, the correct spectrum must be below 
Pao’s form for K 0 = 2.45. The skewness would not change in this case 
because the major contribution to it comes from the wavenumber regime al- 
ready simulated. Experimental spectra by Champagne et al. (1977) show 
a high wavenumber bump which is similiar to the computed spectral shape. 
Additional experiments which determine both the Kolmogorov constant and 
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the skewness and make detailed comparisons with Pao’s spectral form would 
help answer these questions. 

There are no experiments or phenomenological models with which to com- 
pare our results on the dependence of the mixed-derivative skewness with 
Reynolds number. But there are arguments that the mixed-derivative skew- 
ness should increase with Prandtl number, possibly with an upper bound at 
high Prandtl number for given Reynolds number (Hill 1978). Hill’s arguments 
assume that there are two spectral regimes for the scalar variance at low 
Prandtl numbers, a k^ inertial regime and a k^ 1 inertial-diffusive regime. 
There are strong theoretical arguments in favor of both spectral regimes 
(Corrsin 1951 and Batchelor et al. 1959). Experimentally, a k ^ inertial 
range for the scalar variance has been well demonstrated (Monin and Yaglom 
1975, p.511), but the evidence for a A: - ^" inertial-diffusive regime is tenuous. 
This is due in part to the exotic nature of low-Prandtl-number materials, 
such as liquid mercury and sodium. Experimental results on the dependence 
of the mixed-derivative skewness on Prandtl number are also difficult to ob- 
tain. With the exception of a few measurements in air ( Pr — 0.7 in figure 4) 
the only measurements of the mixed-derivative skewness were taken in water 
and liquid mercury by Clay (1973). For air, measurements in the atmosphere 
(Antonia and Van Atta, 1978) give much larger values than those in wind 
tunnels. Our results, Hill’s model #4, the test field model (Larcheveque et 
al. 1980), and several experimental values (taken from Larcheveque et al.) 
are shown in figure 5. Our values are the maximum S u $ for a given Prandtl 
number in our simulations (see figure 3). The experiments of Clay (1973) 
suggest that there is virtually no dependence on Prandtl number, which is 
compatible with our results. Figure 3 suggests that for higher Reynolds num- 
bers S u e will asymptote to -0.5 for Pr = 0.1, but remain constant for Pr = 
0.5 and 1.0. In this case, simulations at higher Reynolds numbers should show 
that the mixed-derivative skewness is constant over a wide range of Prandtl 

g 

numbers. A long crossover regime between the k 3 and k 3 spectral ranges 
could explain this result. Gibson (1968) has suggested that this new regime 
should go as A: -3 , but the spectrum for Pr=0.1 in figure 2 does not support 
this. 

§E. Fourth- and Higher-Order Velocity Correlations 

Siggia (1981b) has shown that in isotropic flow all fourth-order velocity- 
derivative correlations can be expressed in terms of four rotational invariants 
of the velocity-deformation tensor. These are 
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(29) 
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and 

/ 4 = < a/ 4 > 


(31) 


where is the rate of strain, e 2 

, Ui = V x Ui is the vorticity, and w 2 = E w * 2 - 
correlations as follows: 
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and 
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(35) 


Their uncorrelated values are 3, 3, 1, and 3, respectively, and they are plotted 
in figure 6. 

The only fourth-order velocity-derivative correlation that has been ex- 
perimentally determined is the velocity-derivative flatness, Fq u * (3), which 
equals F\ if isotropy is assumed. Kuo and Corrsin (1971) find that Fgu 4 in- 
creases with Reynolds number as R \ a * with a 4 = 0.25 for 12 < R\ < 100. 
Our flatness (Fi or Fq u 4 , figures 6, 7 and 8) also increases with Reynolds 
number and is in excellent agreement with the experiment of Frenkiel and 
Klebanoff (1975) (upper dotted line, figure 7), but a 4 = 0.18 ± 0.03. 

If only one scaling exponent is necessary to describe intermittency (as 
suggested by the phenomenological models), all the fourth-order velocity- 
derivative correlations should exhibit similiar behavior. Correlations Fi to 
F 4 have been calculated for the forced runs and plotted as functions of 
Reynolds number in figure 6. If power laws are fitted to F\ , F 2 , and F 4 
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for R\ greater than 28, the exponents are 0.18 ± 0.03 , 0.29 i 0.03, and 
0.37 i 0.03 , respectively. This suggests that not one, but at least two 
scaling exponents are necessary, one for the rate of strain and the other for 
the vorticity. Furthermore, F 3 decreases noticeably with Reynolds number at 
very low Reynolds number and rises only a small amount at higher Reynolds 
numbers. The low value for F 3 suggests that the principal component of the 
rate of strain and the vorticity tend to be perpendicular to each other, which 
would always be the case in two dimensions, and implies that there is strong 
alignment in the small-scale turbulent structures. This is demonstrated by the 
graphics discussed in section G. The similiarity between one of the isotropic 
forms of the skewness (17b) and / 3 suggests that the low value for F 3 is related 
to the invariance of the velocity-derivative skewness. 

Although the computed fourth-order correlations suggest that several scal- 
ing exponents are necessary, it is possible that at much higher Reynolds num- 
bers one scaling exponent is sufficient. If so, the self-similiar regime described 
by the models has not been reached. The slight rise in F 3 might indicate 
that the calculations are approaching this regime. A change in the velocity- 
derivative flatness scaling exponent at higher Reynolds numbers would be 
seen as all the exponents approach the asymtotic value. This could be in- 
terpreted as a variation in the dissipation-dissipation correlation exponent // 
with Reynolds number. Experimentally, the Reynolds-number dependence of 
the velocity-derivative skewness and flatness is different at very high Reynolds 
number, with a 4 ~ 0.32 for R\ = 200 to 20,000 (Van Atta and Antonia 
1980). There are two experimental measurements that would help clarify 
this matter. If U 2 and u\ can be measured simultaneously, which now seems 
possible with crossed-wire probes, then 



could be measured. If the trends calculated are correct, 1% should dominate 
in I 5 , and / 4 should dominate in Iq. This implies that the scaling exponent of 
the normalized correlation of / 5 should be greater than that of Fg u 4, and the 
scaling exponent of the normalized correlation of should be the largest. 

The fifth- through eighth-order velocity-derivative correlations (3) and vor- 
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ticity moments, 


< ui n > 


(38) 


F^n 


< a ;/ 2 > 


n/2 


have also been calculated. Since the statistics deteriorate as the order of 
the correlations increases (appendix A), only the velocity-derivative skewness, 
flatness, and fifth- and sixth-order correlations are plotted in figure 7, even 
though the fifth- and sixth-order statistics are questionable. The exponents 
for the velocity-derivative correlations in figure 7 are 0.18 d: 0.03, 0.27 ± 
0.03, and 0.55 ± 0.03 respectively. The lognormal theory predicts that 
these exponents should be proportional to n(n-l) (equation 4) and the /?- 
model predicts (n-2) (equation 5). The lognormal theory provides a better 
fit, as in the experiments (Antonia et al. 1982). The associated value of 
the dissipation-dissipation correlation exponent // has not been calculated 
because the simulations are not in the Reynolds number regime where the 
experimental value is usually found. 

Assuming isotropy, F u 4 = F4. The most important point that can be 
made about the sixth-order vorticity correlation (F w e, figure 9) is that its 
scaling exponent (1.1 i 0.1) is much larger than the corresponding sixth-order 
velocity-derivative scaling exponent, just as was true for the fourth-order 
correlations. 


§F. Scalar and Mixed Scalar- Velocity Fourth-Order Derivative Correlations 

There are three fourth-order derivative correlations which have been 
measured experimentally, the velocity-derivative flatness, Fq u 4 (3), the scalar- 
derivative flatness 


Fd»n = 


< (&r > 

< (rer ) 2 >" /2 ' 


n = 4 


and the mixed-derivative correlation 


(39) 


F u g = 


<< ( dj \2 
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(40) 


The uncorrelated value for the velocity- and scalar-derivative flatnesses, as- 
suming Gaussian statistics, is 3, and the uncorrelated value for the mixed- 
derivative correlation is 1. The rotationally-invariant form of the scalar- 
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derivative flatness in isotropic flow is 


Fd $ 4 — Fve = 


9 < V0 4 > 

5 < ve 2 > 2 


(41) 


In the experiment of Antonia and Chambers (1980), the atmospheric 
temperature-derivative flatness increases much faster with Reynolds num- 
ber than the velocity-derivative flatness. Specifically, if Fqo* ~ R\ a * , then 
ote = 0.5 for 100 < R\ < 10,000. In our simulations (figure 8), the scalar- 
derivative flatness also has a stronger dependence on Reynolds number than 
the velocity-derivative flatness. In addition, the scalar-derivative flatness in- 
creases with Prandtl number for given Reynolds number, while the scaling 
exponents decrease. For Prandtl numbers 0.1, 0.5, and 1.0 we found ae = 
0.48 , 0.44, and 0.36. The trend for the flatness factors at low Prandtl num- 
bers to increase faster with Reynolds number than at larger Prandtl numbers 
suggests that the large exponents for Pr = 0.1 and 0.5 might be transient, 
and that at higher Reynolds numbers they will have a exponent similiar to 
that for Pr — 1.0. That is, a$ = 0.36, which is the same order as the ex- 
ponent for the vorticity flatness (P 4 ), might be an upper bound for the scaling 
exponent at large Reynolds numbers for all Prandtl numbers. 

Although the scalar-derivative flatness increases with R\ at similiar rates in 
experiment and in our simulations, the Reynolds numbers are very different. 
Sreenivasan et al. (1980) have measured the scalar-derivative flatness in a 
wind tunnel with temperature fluctuations produced by a heated screen and 
find Fqo * = 5.5. The Prandtl number is 0.7 and we believe R\ m 24. (fix 
was computed from their u 2 data and is consistent with similiar wind-tunnel 
experiments by Warhaft and Lumley (1978).) For R\ & 24, their result is 
consistent with our calculations and lies on our curve for Pr = 0.5 in figure 
8 . 

The sixth-order sc alar- derivative correlations, Foe 6 , are plotted in figure 9. 
The scaling exponents for Prandtl numbers 0.1, 0.5, and 1.0 are 1.48, 1.30, 
and 0.92 respectively, showing the same trends with Prandtl number as the 
fourth-order correlations. 

Before disscussing the one experimentally-measurable mixed correlation we 
will consider some rotationally-invariant correlations between the velocity and 
scalar derivatives where both are second order. These are 


_ < (vo)V > 

V * wl <V9 2 Xu 2 > 


(42) 


20 


(43) 


Fv$u2 = 3 


< (Vfl • o >) 2 > 
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Fv del = 
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and 

3 < de/dxieije jk d9/dx k > 

Fv,e2 = < V0 2 > < e 2 > (45) 

They have been normalized so that their uncorrelated values are 1. 
Correlation Fv^wi is plotted in figure 10 and is nearly 1, showing that there is 
little correlation between the magnitudes of the vorticity and scalar deriva- 
tive. Correlation Fveu 2 is plotted in figure 11 and shows a strong anti- 
correlation between the direction of vorticity and the direction of the scalar 
gradient, which would be consistent with the scalar is being wound around 
vortex tubes. Batchelor (1952) predicts this result for the stretching of a 
volume element of fluid. 

Both correlations between the strain rate and the scalar derivative are 
greater than one and show a tendency to increase with Reynolds number for 
R\ < 30, then remain constant or increase slowly, much like the skewnesses 
do. The experimentally-measurable mixed correlation is related to F^e el and 
Fvoe 2 by 


< (|2Ll)2(^4_)2 >0 4 

F “' = <(%?><, ftp > = + 7*- 

and is plotted in figure 12(a). F u $ has been measured by Park (1976) and is 
greater than one. Correlations F^ee 1 and 2 are plotted in figures 12(b) 
and 12(c) respectively. 

We have calculated a wide assortment of other correlation factors, none of 
which show a significant divergence from their uncorrelated values. This is 
especially true for those with the scalar variance, such as 

< 9 * > 

< 9 2 > 2 

For both our simulation and experiment (Sreenivasan et al. 1980) the scalar- 
variance flatness is 3, the Gaussian value. This is expected because, except for 
mixing, the scalar level is not affected by the dynamics although its derivatives 
are. 
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§G. Graphics In Siggia (1981a) three-dimensional perspective plots of the 
vorticity field indicated the presence of vortex tubes. We have produced 
similiar graphics for the vorticity, scalar gradient, and principal rates of 
strain for a single time from run F19. These are shown in figures 13 to 18. 
These figures are two-dimensional, black-and-white projections of what can 
be seen using colored graphics or a three-dimensional graphics terminal. In 
figures 13(a-c) short lines are used to show the direction of vorticity, with 
the length of the lines proportional to the magnitude of the vorticity. The 
vorticity is plotted at a grid point only if it was above a threshold, which is 
picked such that only 0.25%, or about 500, of the grid points of a 64 3 mesh 
are plotted. The 0.25% value was picked so as to allow an assortment of 
structures to appear, but not so many that the pictures would be cluttered. 
Figure 13(b) represents the same flow field as figure 13(a), but rotated 90° 
about the vertical axis. A strong concentration of vorticity extends most of 
the way across our computational box in both projections. Rotation shows 
that this structure is a tube: The vorticity in the gap indicated by the circle 
in figure 13(a) is just below the threshold picked. In figure 13(c) the threshold 
is lowered such that 1% of the grid points are plotted and this gap is filled. 
The circled area in figure 13(c) indicates where two tubes might be merging 
into one. 

In figures 14(a,b) the lines represent the direction of the scalar gradient. 
The projections are the same as in figures 13(a,b) and 0.25% of the grid points 
are plotted. The scalar gradient field appears to have well-defined structures, 
and although it is not very clear in figures 14(a,b), rotation shows that these 
structures are sheets. By comparing figures 13(a) and 14(a) we can see that in 
the vicinity of the largest vortex tube that the scalar gradient is perpendicular 
to the vorticity. This is consistent with our calculation of 2 (43), which 
suggested that the scalar gradient and vorticity are orthogonal. 

To better visualize the correlation between vorticity and scalar gradient, 
two-dimensional slices of the vorticity and scalar gradient are plotted in 
figures 15 and 16. The location of the slices is indicated by the arrow in figure 
13(a). In these pictures, the scalar gradient (figure 16) appears wound around 
the vortex core (figure 15). Two-dimensional vortex simulations which include 
scalar markers (Aref and Siggia 1980) and flow-visualization experiments 
(Winant and Browand 1974; Brown and Roshko 1974) also show the scalar 
marker winding around what are believed to be vortex cores. 

In figures 17(a,b) the rate of strain is plotted for the same flow field and 
projections as in figures 13(a,b). Since the rate of strain is a tensor, it cannot 
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be described with a single line segment. However, by plotting only one of 
the principal rates of strain at a time, graphics can be produced that show 
the structure of the strain field. In figures 17(a,b) we conditionally sample 
based on the trace of the square of the rate-of-strain tensor (e 2 ), then plot 
the principal rate of strain with the largest absolute value; 0.5% of the grid 
points are used and in every case the largest principal rate of strain was 
negative, or compressive, with stretching in the two directions perpendicular 
to this component to maintain incompressibility. We might expect the largest 
principal rate of strain to be compressive because the velocity-derivative 
skewness (14) is negative. That is, by comparing the skewness with one of 
its isotropic forms (the trace of the rate of strain cubed, (17a)), we see that 
a negative skewness implies a negative value for the largest principal rate of 
strain. 

In figures 17(a,b) there are obvious structures in the rate-of-strain field, but 
it is not clear from rotation whether these are tubes or sheets. However, there 
is definite alignment with respect to the vorticity and scalar gradient which is 
consistent with the calculated statistics. First, because F% (33), the correla- 
tion between the magnitudes of the rate of strain and vorticity, is large, we 
expect the rate of strain to be concentrated near the vortex tubes. In figure 
17(a) the largest strain structure is located near the strong vortex tube in 
figure 13(a). Next, since there will, on the average, be vortex stretching along 
the tubes, any compressive components should be perpendicular to the tubes. 
By comparing figures 13(a) and 17(a) it appears that the compressive com- 
ponent of the rate of strain is aligned perpendicular to the vortex tube. This 
is further illustrated by comparing figure 18, which plots a two-dimensional 
slice of the strain, with figure 15, the slice of the vorticity. The slice in figure 
18 is the same as in figures 15 and 16. Finally, how do the stretching, or ex- 
pansive, components of the rate of strain align with the vorticity? A linearized 
model due to Viellefosse (1982) shows the larger expansive component aligned 
perpendicular to the vorticity, and the smaller expansive component aligned 
along the vorticity. Since one of the rotationally-invariant forms of the skew- 
ness (17b) is related to vortex stretching, and both the skewness and F s are 
small, our statistics also suggest that the smaller expansive principal rate of 
strain is aligned with the vorticity. Unfortunately, we have no graphics that 
show the expansive components of the rate of strain clearly. 

How would we expect the scalar gradient to be aligned with respect to 
the rate of strain? First, the scalar-gradient sheets and rate-of-strain struc- 
tures should be located near one another because the fourth-order correla- 
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tions between the rate of strain and scalar gradient (44,45) are greater than 
one. Second, the scalar gradient should be aligned along the compressive 
component of the rate of strain because the mixed-derivative skewness is 
negative. That is, just as we predicted that the largest component of the rate 
of strain is compressive by comparing the velocity-derivative skewness with 
one of its isotropic forms, by comparing the mixed-derivative skewness with 
its isotropic form (18) we predict that the compressive component of the rate 
of strain is aligned with the scalar gradient. Comparisons between figures 
14a, b and 16 to 17a, b and 18 show that both of these predictions are borne 
out. 

§H. Discussion We have studied the derivative correlations of an idealized 
turbulent flow because they are a quantitative means of probing the small 
scales of that flow. Numerically, the derivative correlations helped establish 
the resolution of our simulations. Once the resolution was established, they 
were compared with experiments and phenomenological models and were used 
in the interpretation of three-dimensional graphics. Because the small scales 
of turbulence are believed to be decoupled from the large scales, what is 
learned about the small scales in this idealized flow should be applicable to 
turbulent flows in general. 

The most important result for the velocity-derivative correlations is that 
they do not obey the scaling laws predicted by statistical models of intermit- 
tency, such as the lognormal model and the /?-model. These models predict 
that the skewness will depend on Reynolds number and that all the fourth- 
order velocity-derivative correlations will have the same Reynolds-number 
dependence. In the simulations the skewness does not depend on the Reynolds 
number and each of the rotationally-invariant fourth-order correlations has 
a different exponent. It is possible that at very high Reynolds numbers the 
velocity-derivative correlations do behave as these models predict and there 
are high-Reynolds- number experiments (Tavoularis et al. 1978) that show 
the skewness increasing and a different power law for the velocity-derivative 
flatness (Van Atta and Antonia 1980). Alternatively, models that include the 
structure of the small scales, such as the model of Tennekes (1968), might be 
necessary. The skewness is consistent with the suggestion of Tennekes that 
vortex stretching is caused by large eddies the size of the Taylor microscale. 
Tennekes’s model also suggests that the velocity-derivative flatness increases 
with Reynolds number with <24 = 1.0. Both our simulation and experiments 
find a much smaller value for a 4 . In appendix B it is shown how simple 
modifications of Tennekes’s model can give smaller values of a 4 , while 0:3 
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remains zero. This might provide a basis for a more complete model of small- 
scale turbulence. A simple experiment which might resolve this issue is to 
find the exponents of more of the fourth-order correlations. 

The most important result for the derivative correlations of a passive scalar 
is that the scalar-derivative flatness is larger than the velocity-derivative 
flatness. This suggests that the scalar derivative is more intermittent than 
the velocity derivative and might be related to several other anomalous scalar 
effects seen experimentally, for example the “bump” in the scalar spectrum 
described by Hill (1978) and sharp gradients, or interfaces, called “ramps” 
(Antonia et al. 1979). What allows sharp structures to form in the scalar 
field, but not the velocity field? First, they do appear in the velocity field, but 
as vorticity, which cannot be measured directly, and not as simple velocity 
gradients, which are easy to measure. These structures are much more likely 
to be observed in the du 2 fdx i field, than in duifdx\, if our conclusions about 
the relative strengths of 1 5 and Iq (36,37) are correct. 

Graphical display has been used to demonstrate the alignment indicated 
by the statistics. This shows that the vorticity is concentrated in tubes, not 
sheets, with large concentrations of the rate of strain and scalar gradient 
nearby. The statistics and graphics show that the largest principal rate 
of strain is compressive and aligned perpendicular to the tube. The larger 
stretching, or expansive, component of the rate of strain also appears to be 
perpendicular to the vortex tubes, while the stretching along the tubes is small 
and is probably not caused by the immediate vortex tube. This structure and 
the statistics are consistent with the linearized model of Viellefosse (1982) 
and the structure model of Tennekes (1968), who predicted that the velocity- 
derivative skewness is independent of Reynolds number. Graphical display 
also shows that large values of the scalar gradient are wound in sheets around 
the tube and that the gradient is aligned perpendicular to the vorticity and 
along the compressive component of the rate of strain. This is consistent with 
correlations between the scalar gradient and both the rate of strain and the 
vorticity. 

Although flow visualization has shown the existence of extended “coherent” 
structures on the large scales, most theoretical models of the small scales of 
turbulence have neglected alignment. They assume that small-scale struc- 
tures in fully-developed turbulence are distributed with random orientation 
and that statistical mechanics is applicable. Our results suggest that tur- 
bulence is characterized by extended vortex tubes and strong alignment bet- 
ween the vorticity and rate of strain. It is possible that as the Reynolds 
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number increases these tubes become more tightly wound and the alignment 
becomes stronger. There would still be statistics, but the statistics would 
have to include structures. Strong alignment is not inconsistent with the best 
known characteristic of turbulence, the A; -5 / 3 spectrum. Lundgren (1983) has 
shown how fluctuations about a Burger’s vortex can give a A; -5 / 3 spectrum. 
Our highest Reynolds number simulation also shows a short -5/3 regime. 
Even if the correlations behave as the statistical phenomenological models 
predict at higher Reynolds numbers, there should still be a tendency for the 
vorticity to form tubes and for there to be alignment between the vorticity 
and rate of strain. The vortex tubes might show more random orientation 
than in our current graphics, but there would still be strong alignment at the 
smallest scales. 

Often, subgrid modelling is used with a three-dimensional code (Moin and 
Kim 1982) to achieve higher Reynolds numbers. This might be useful for 
studying scalar transport, but would be inappropriate in our work. Subgrid 
modelling would smooth the flow and result in lower, and incorrect, skewness 
and flatness factors. 

The next step in studying small-scale statistics with numerical codes will 
be the influence of large-scale shear, strain and temperature gradients. For 
example, it is known that in the presence of a shear and temperature gradient 
that the sc alar- derivative skewness (16) is nonzero. We also want to learn 
more about the universality of the isotropic statistics in anisotropic flows and 
study “band- averaged” statistics. By band-averaged we mean calculating 
flatnesses based on only part of the velocity spectrum. This has been done 
experimentally (Tavoularis and Corrsin 1981) and comparisons will be made. 
After that we will experiment with different subgrid models and see how well 
they preserve the band-averaged statistics. We hope that once a clear picture 
of the small-scale structures is developed that improved methods of subgrid 
modelling can be found. 

Future areas of work with passive scalars are the return to isotropy of 
an anisotropic scalar distribution, the dispersion of an interface, the correla- 
tion between two independent scalars, and, eventually, chemical reactions. 
Turbulent chemical reactions are believed to occur primarily along the sharp 
interfaces that develop in the scalar field. Intermittency will have a significant 
effect on the development of these interfaces. 
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§Appendix A: Errors Many of the trends calculated in this paper, in particular 
the dependencies on Reynolds and Prandtl numbers, could be sensitive to 
moderate to large errors in our statistical sample. Therefore, in tables 3-6 
the correlations and their statistical variances are given. Where possible, 
isotropic forms of the correlations have been used to effectively triple the 
sample size and reduce the errors. For one case, F7b, errors are given based 
on comparisons with a much longer simulation. By comparing the variances 
and errors in F7b we can determine a characteristic sample size for each 
statistic, N s , which we define by 

/ v (VQrTiCL7lCC)§ 

y*pCT C&7ltd£l £ 6TT OT ) 6 — ■ . 

To calculate the errors for the other forced simulations, first estimate the 
sample size by comparing the number of eddy-turnover times with the number 
in F7b. Then calculate the error by 

. . (variance) s 

(percentage error) s = — , 

n/aW/aW - 1 


where Np 7 is the number of eddy-turnover times in F7b, N e is the number 
in F13, and N s is the sample size for the flatness in F7b. As an example, we 
find the error for the velocity-derivative flatness for F13 (table 3), 


2 . 0 % = 


5.2% 

x /l4(2.6/4.5) - 1 


Detailed tables of the errors for the 32 3 simulations are given in Kerr (1981). 
All the third- and fourth-order correlations have error bars which are less than 
10%. Fifth- and higher-order correlations and variances were not calculated 
for F7, but based on the large variances for these statistics in the other 
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simulations (tables 5 and 6), we expect that their errors are very large. Fifth- 
and higher-order correlations are included only to allow a rough comparison 
with the phenomenological models. Third- and fourth-order statistics in 
tables 5 and 6 were calculated without the isotropic forms discussed and as 
a result their errors are much higher than those shown in tables 3 and 4. 

Because of their significance, the time-dependence of the velocity-derivative 
skewness (14) and one of the fourth-order velocity-derivative correlations, 
F3 (34), are shown for several Reynolds numbers in figures 3(a) and 6(a), 
respectively. The time-scale has been normalized by the large-scale eddy- 
turnover time, t e , for each Reynolds number. These figures demonstrate that 
the trends discussed in two earlier sections (Spectra and Skewness; Fourth- 
and Higher Order Velocity Correlations) are not statistical errors. Typical 
times over which samples were taken are indicated by the dashed lines in 
figures 3(b) and 6(a). 

§Appendix B: Extension of Tennekes model Let us first review the assump- 
tions made by Tennekes (1968) to find the Reynolds-number dependence of 
the velocity-derivative skewness and flatness. His first assumption was that 
the small-scale vortices are strained by vortices whose scale is the Taylor 
microscale, X. This is equivalent to assuming that the small-scale vortices 
feel only the average rate of strain, 


Next, he assumed that the small-scale vortices have a characteristic width, 
which we will call f. He assumed that f = r], the Kolmogorov microscale, but 
we will only assume that X > f > rj. This is equivalent to assuming that it is 
some sort of average vortex tube that is significant. To find the higher-order 
correlations, Tennekes assumed that only the vorticity within the small-scale 
structures matters. This is the average vorticity squared, (U/\) 2 , divided by 
the volume occupied by the small-scale structures, which if they are tubes of 
width $ and bend over a distance X, is 


U/\. 


(Bl) 


(f/X ) 2 


(B2) 


Therefore, within the tubes the vorticity is 


(U/Xf = U 

(f/X ) 2 f 


(B3) 
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To find the skewness, Tennekes used (1Tb) and multiplied the vorticity within 
the tubes (B3), the average rate of strain (Bl), and the volume of the vortices 
(B2), then normalized by ( U /X) 3 . This gives 


S u ~ 


< ueu > 
(e/i/) 3 / 2 




(C//f) 2 ([//X)(f/X) 2 

(U/\f 


Notice that all factors of f and X drop out, so it is not critical that the size 
of the small-scale structures be assumed to be the Kolmogorov microscale. 

To find the flatness, one simply squares the vorticity within the small-scale 
structures, multiplies by their volume, and normalizes by (f7/X) 4 . This gives 


Fdu 4 — 


(r//c ) 4 (f/x ) 2 

(U/\f) 



If g = n, then as Tennekes showed, Fg u 4 ~ R\, or a 4 = 1. In our 
simulations, with a 4 < < 1, we require that 

f = X x Rx aA/2 


It also follows that when n is even that ot n = q 4 x rk ^. This is the 
same equation, (5), as found by the ^-model of Frisch, Sulem, and Nelkin 
(1978). When n is odd our model suggests that a n = a( n _i). As pointed 
out, our simulations and experiments do not support the /?-model. Nor does 
this simple model explain why the velocity-derivative flatness and vorticity 
flatness have different scaling exponents. But a simple model such as this 
might be the starting point of a more complete model of the small scales. 
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TABLE 1.- CHARACTERISTICS OF FORCED SIMULATIONS 


Fun 

Mesh 

r A 

Pr 

At a 

N f fc 

*k 

KOC 

t e 

FI 

32 3 

9.0 

2.0 

19.0 

20 

9.4 

15.8 

4.7 

F2 

32 3 

9.0 

. 5 

19.0 

20 

9.4 

5.6 

4.7 

F3 

32 3 

12.3 

2.0 

13.6 

35 

11.9 

20.0 

3.55 

F4 

32 3 

12.3 

1.0 

13.6 

35 

11.9 

11.9 

3.55 

F5 

32 3 

12.3 

. 5 

13.6 

35 

11.9 

7.8 

3.55 

F6 

32 3 

18.5 

1.0 

9.6 

25 

16.5 

16.5 

1.97 

F7a 

32 3 

18.5 

. 5 

39.0 

40 

16.4 

9.8 

1.98 

F7b 

32 3 

18.5 

.5 

9.0 

10 

16.5 

9.8 

1.97 

F8 

32 3 

18.5 

. 1 

9.6 

25 

16.5 

2.9 

1.97 

F9 

32 3 

24.0 

. 5 

9.5 

20 

21.4 

12.7 

1.26 

F10 

32 3 

24.0 

. 1 

9.5 

20 

21.4 

3.8 

1.26 

FI 1 

32 3 

28.9 

. 5 

29.5 

60 

25.7 

15.3 

. 91 

Fl 2 

32 3 

28.9 

— 

29.5 

60 

25.8 

— 

. 91 

F13 

64 3 

28.5 

1.0 

10.0 

11 

12.3 

12.3 

3.83 

F14 

64 3 

28.5 

. 5 

10.0 

11 

12.3 

7.3 

3.83 

F15 

64 3 

28.5 

. 1 

10.0 

11 

12.3 

2.2 

3.83 

F16 

64 3 

37.5 

1.0 

4.5 

10 

16.0 

16.0 

2.51 

F17 

64 3 

37.5 

. 5 

4.5 

10 

16.0 

9.5 

2.51 

F18 

64 3 

37.5 

. 1 

4.5 

10 

16.0 

2.8 

2.51 

F19 

64 3 

48.2 

1.0 

6.75 

28 

22.4 

22.4 

1.41 

F20 

64 3 

48.2 

. 5 

6.75 

28 

22.4 

13.3 

1.41 

F21 

64 3 

48.2 

.1 

6.75 

28 

22.4 

4.0 

1.41 

F22 

64 3 

55.9 

1.0 

4.00 

17 

27.2 

27.2 

1.05 

F23 

64 3 

55.9 

. 5 

4.00 

17 

27.2 

16.2 

1.05 

F24 

64 3 

55.9 

. 1 

4.00 

17 

27.2 

4.8 

1.05 

F25 

128 3 

82.9 

1.0 

1.35 

27 

45.7 

45.7 

.85 

F26 

128 3 

82.9 

. 5 

1.35 

27 

45.7 

27.2 

. 85 

F27 

128 3 

82.9 

. 1 

1.35 

27 

45.7 

8.1 

. 85 


^Time span of statistical sample. 
Number of files averaged. 
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TABLE 2 


VELOCITY AND SCALAR FLATNESS FACTORS USED TO 
DETERMINE ALIASING ERRORS 


Run 

Grid size 

Pr 

F 3u 4 (6 4 ) 

F 3e ,(6 4) 

F 3u 4(32) 

F 30 4 < 32 > 

Averages of many realizations 

Fll 

32 

0.5 

— 

— 

3.94 

5.73 

Fl 2 

32 

— 

— 

— 

4.17 

-- 

F13 

64 

1.0 

4.18 

6.88 

-- 

-- 

F14 

64 

.5 

4.18 

5.82 

-- 

— 

F15 

64 

. 1 

4.18 

3.61 

— 

— 

One realization 

F13 

64 

1.0 

2.99 

7.30 

3.93 

5.62 

F14 

64 

. 5 

3.99 

5.95 

3.93 

5.08 

F15 

64 

.1 

3.99 

3.55 

3.93 

3.57 


Note: FI 1 is 32 3 aliased. F12 is 32 3 de-aliased. The size 

of the grid used to calculate the flatness is in parentheses. 


TABLE 3.- VELOCITY DERIVATIVE SKEWNESS AND FLATNESS FACTORS 

AND VARIANCES 


Run 

r a 

S u 

s 

e 

F i 

F 4 

F 2 

F 3 

Decaying 

2.33 

-0.32 






Decaying 

4.07 

-.40 






Decaying 

6.20 

-.44 






F1-F2 

9.0 

-.502 

0.868 

3.39 

3.17 

2. 26 

0.693 



6.7% 

1.6% 

5.3% 

4.9% 

6.5% 

5.3% 

F3-F5 

12.3 

-.523 

.748 

3.51 

3.35 

2.51 

. 642 



6.9% 

1.5% 

4.6% 

6.5% 

6.5% 

4.6% 

F7a 

18.5 

-.506 

. 621 

3.69 

3.82 

3.72 

. 612 



3.8% 

. 2% 

5.4% 

7.9% 

7.7% 

7.7% 

F7b Real error 


.7% 


1.5% 

1.8% 

1.7% 

. 6% 

F9-F10 

24.0 

-.507 

. 578 

3.92 

4.36 

4.14 

.621 



2.2% 


3.7% 

8.1% 

6.0% 

5.9% 

F13-F15 

28.5 

-.507 

.553 

4.18 

4.53 

4.40 

. 658 



.0% 


5.2% 

8.3% 

7.7% 

5.7% 

F16-F18 

37.2 

-.504 

.530 

4.26 

4.86 

4.53 

. 652 



2.4% 


4.3% 

7.8% 

5.6% 

6.4% 

F19-F21 

48.2 

-.501 

. 521 

4.38 

5.42 

4.88 

.675 



1.7% 


3.5% 

5.2% 

4.9% 

4.5% 

F22-F24 

55.9 

-.497 

.520 

4.54 

5.77 

5.08 

.712 



1.6% 


4.4% 

6.1% 

4.9% 

4.9% 

F25-F27 

82.9 

-.511 

.523 

5.11 

7.09 

6.00 

.779 



. 9% 


3.8% 

5.9% 

4.6% 

3.6% 
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TABLE 4.- SCALAR- AND MIXED -DERIVATIVE SKEWNESS AND 
FLATNESS FACTORS AND VARIANCES 


Run 

r a 

Pr 

S u6 

s x 

<1 

CD 

F V0o)i 

F 

V0ei 

F 

V0OJ2 

F 

V0e 2 

FI 

9.0 

2.0 

-0.523 

3.8% 

0.585 

5.18 

14.3% 

0.94 

5.1% 

1.16 

6.1% 

0.663 

10.4% 

1.77 
7. 9% 

F2 

9.0 

.5 

-.420 

3.3% 

. 749 

3.74 

5.8% 

. 97 
4.6% 

1.07 

4.4% 

.768 

11.7% 

1.49 

6.4% 

F3 

12.3 

2.0 

-.520 

3.8% 

.554 

5.83 

11.4% 

. 93 
5.6% 

1.18 

6.4% 

.636 

8.8% 

1.82 

8.0% 

F4 

12.3 

1.0 

-.510 

3.9% 

. 608 

5.11 

11.2% 

. 93 
6.5% 

1.16 

6.8% 

.657 

8.0% 

1.76 

8.5% 

F5 

12.3 

.5 

-.465 

4.3% 

. 673 

4.25 

6.9% 

. 95 
5.7% 

1.12 

5.8% 

.708 

8.3% 

1.64 

7.5% 

F6 

18.5 

1.0 

-.506 

4.9% 

. 551 

5.87 

9.4% 

. 92 
5.3% 

1.19 

6.0% 

. 624 
7.5% 

1.82 

7.1% 

F7 a 

18.5 

.5 

-.479 

3.1% 

.579 

4.93 

7.2% 

. 94 
8.3% 

1.12 

4.0% 

.690 

8.8% 

1.69 

4.0% 

F7b 

Real errors 

±. 6% 


-2.7% 

+ 2.0% 

+ 2.8% 

±. 9% 

+ 2.8% 

F8 

18.5 

.1 

-.286 

8.7% 

. 813 

3.37 

2.9% 

. 98 
4.9% 

1.03 

4.9% 

. 843 
5.4% 

1.31 

5.3% 

F9 

24.0 

.5 

-.486 

4.1% 

. 549 

5.63 

12.2% 

. 94 
5.1% 

1.18 

5.3% 

. 663 
7.9% 

1.81 

5.8% 

F10 

24.0 

.1 

-.344 

7.3% 

. 657 

3.58 

6.0% 

. 98 
3.9% 

1.05 

4.0% 

. 801 
5.6% 

1.43 

4.6% 

F13 

28.5 

1.0 

-.517 

3.6% 

.527 

6.88 

8.4% 

1.03 

7.2% 

1.38 

5.8% 

.217 

9.6% 

2.14 

6.2% 

F14 

28.5 

. 5 

-.506 

3.9% 

.544 

5.82 

10.9% 

1.02 

6.9% 

1.33 

4.6% 

. 281 
9.9% 

2.05 

4.9% 

F15 

28.5 

. 1 

-.361 

5.7% 

.596 

3.61 

10.2% 

1.00 

7.3% 

1.13 

4.2% 

.557 

6.7% 

1.56 
4 . 6% 

F16 

37.2 

1.0 

-.521 

2.5% 

.523 

7.57 

7.5% 

1.06 

6.6% 

1.41 

5.6% 

. 154 
8.9% 

2.21 

5.7% 

F17 

37.2 

.5 

-.517 

2.7% 

.535 

6.49 

6.3% 

1.06 

6.1% 

1.39 
4 . 9% 

. 205 
8.8% 

2. 17 
4.9% 

Fl 8 

37.2 

.1 

-.408 

3.7% 

.549 

3.90 

7.1% 

1.06 

5.0% 

1.20 

4.0% 

.464 

6.6% 

1.73 

3.8% 

FI 9 
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.5 
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.532 
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10.7% 
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4.2% 
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4.6% 

F21 

48.5 

.1 

-.446 

4.3% 
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TABLE 5.- HIGHER-ORDER VELOCITY DERIVATIVE AND VORTICITY CORRELATIONS 

AND VARIANCES 


Run 

R X 

S 

u 

F 3u 4 

F 3u 5 

F 3u 6 

Fr> 7 

dU 

F 3u 8 

F ^ 

0 J 

V 

F 8 

0) 

F13-15 

28.5 

0.527 

4.12 

-6.86 

39.5 

-118 

701 

4.69 

51.8 

988 



24% 

11% 

40% 

41% 

82% 

93% 

13.4% 

40% 

70% 

F16-18 

37.2 

.531 

4.22 

-6.98 

40.1 

-117 

638 

5.00 

61.2 

1302 



21% 

7.8% 

28% 

24% 

45% 

55% 

9.1% 

24% 

39% 

F19-21 

48.5 

.525 

4.38 

-7.61 

47.2 

-159 

1003 

5.49 

82.3 

2620 



20% 

7.7% 

32% 

30% 

83% 

•- 76% 

7.7% 

35% 

88% 

F22-24 

55.9 

.528 

4.63 

-8.07 

53.8 

-179 

1206 

5.85 

96.7 

3284 



17% 

7.4% 

25% 

25% 

50% 

62% 

8.7% 

33% 

76% 

F25-27 

82.9 

.526 

5.11 

-9.38 

73.3 

-281 

2326 

7.12 

185.0 

15613 



10% 

5.2% 

20% 

48% 

47% 

56% 

6.9% 

48% 

142% 


TABLE 6.- HIGHER-ORDER SCALAR- 
DERIVATIVE CORRELATIONS AND 
VARIANCES 


Run 

r a 

Pr 

F ae 4 

F 39 6 

F 8 0 8 

F15 

28.5 

0.1 

3.71 

30 

457 




25% 

67% 

103% 

F18 

37.2 

.1 

3.87 

31 

379 




13% 

33% 

55% 

F21 

48.5 

. 1 

4.50 

46 

784 




22% 

70% 

140% 

F24 

55.9 

. 1 

5.68 

89 

2696 




25% 
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127% 

F27 

82.9 

. 1 

7.32 
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6932 




17% 

48% 

87% 
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28.5 

. 5 

5.87 

89 

2323 




20% 

51% 

82% 

F17 

37.2 

. 5 

6.54 
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3222 
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33% 

56% 

F20 

48.5 

. 5 

7.29 
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4928 




16% 

47% 

84% 

F23 

55.9 

.5 

8.34 
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9300 




18% 

55% 

100% 

F26 

82.9 

.5 

10.8 

371 

24385 




14% 

39% 

66% 

Fl 3 

28.5 

1.0 

6.94 
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4139 




17% 

43% 

72% 

F16 

37.2 

1.0 

7.68 

160 

5669 




15% 

37% 

62% 

F19 

48.5 

1.0 

7.95 

174 

6772 




14% 

41% 

74% 

F22 

55.9 

1.0 

8.02 

183 

7586 




15% 

45% 

85% 

F25 

82.9 

1.0 

10.8 

368 

23398 




13% 

38% 

67% 
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Figure 1.— Kinetic-energy spectra normalized by the Kolmogorov microscales (10a) and multiplied by k 5/3 
Pao’s theoretical curve (13) for K 0 = 2.45 is indicated by the dashed line. Circle: R ^ = 82.9. Triangle 
R x = 55.9. Plus: R x = 37.2. Cross: R x = 1 8.4. 
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Figure 2.— Kinetic-energy and passive-scalar spectra at R ^ = 82.9 normalized by the Kolmogorov microscales 
(10a) and the scalar-variance dissipation, x- Circle: Pr = 0.1. Triangle: Pr = 0.5. Plus: Pr = 1.0. Cross: 
Kinetic energy. 
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(a) Dependence of the velocity-derivative skewness (14) and mixed-derivative skewness (15) on Taylor-microscale Reynolds number. 
Dashed line is the experimental data of Klebanoff (private communication) where S u = 0.43. Square: Pr = 2.0. Circle: Pr = 1.0. Tri- 
angle: Pr = 0.5. Plus: Pr = 0.1. Cross: Kinetic energy. 


Figure 3.— Velocity and mixed-derivative skewnesses. 
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TIME/t e 

(b) Time evolution of the velocity-derivative skewness (14) for three Reynolds numbers. Time has been normalized by the eddy-turnover 
time of each simulation. Square: R ^ = 82.9. Circle: R ^ = 55.9. Triangle: R ^ = 48.2. Dashed line represents the time-span over which 
the ensemble average was taken. 


Figure 3.— Concluded. 
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Figure 4.— Dependence of the velocity-dissipation skewness (24) and the scalar-dissipation skewness (25) on Taylor-microscale Reynolds 

number. 
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Figure 5.— Dependence of the mixed-derivative skewness (15) on Prandtl number. Test field model results (dashed line), the model of Hill 
(1978) (dotted line), and experiments (crossed circles) are all taken from Larcheveque et al. (1980). Symbols for our simulations (solid 
line) are the same as in figure 3a. 
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(a) Dependence on Taylor-microscale Reynolds number. Squares are F, (32) or the velocity-derivative flatness, Fg^4 (3). Circles are F 2 

(33). Triangles are F 3 (34). Pluses are F 4 or the vorticity flatness, F w 4 (38). 

Figure 6.— Rotationally invariant fourth-order velocity derivative correlations. 
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(b) Time evolution of F 3 (34), for three Reynolds numbers. Time has been normalized by the eddy-turnover time of each simulation. 
Square: R ^ = 82.9. Circle: R ^ = 55.9. Triangle: 7?^ = 48.2. Dashed line represents the time-span over which the ensemble average was taken. 


Figure 6.— Concluded. 



Figure 7.- Dependence of the velocity-derivative correlations, Fft u n (3) on Taylor-microscale Reynolds number. Squares are third order, or 
the velocity-derivative skewness, S u (14). Circles are fourth order, the velocity-derivative flatness or F x (32). Triangles are fifth order. 
Pluses are sixth order, which has been divided by 5 so it will appear on this graph. The dotted lines are experimental data of Klebanoff 
(private communication) for the velocity-derivative skewness and the velocity-derivative flatness (stars). 
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Figure 8.— Dependence of the velocity-derivative flatness (3,32) and the scalar-derivative flatness (39.41) on Taylor-microscale Reynolds 
number. The experiment of Sreenvasan et al. 1980 is a star. 
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Figure 9.— Dependence of the sixth-order velocity-derivative correlation, vorticity correlation, and scalar-derivative correlations on Taylor- 
microscale Reynolds number. Squares are the velocity, Fg (/ 6 (3). Circles are the vorticity, f w 6 (38). Triangles are the scalar derivative 
correlation, F$q6 (39), for Pr = 0.1 . Pluses are for Pr = 0.5. Crosses are for Pr = 1 .0. 



Figure 10.— Dependence of the first scalar derivative-vorticity correlation, (42), on Taylor-microscale Reynolds number. 
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Figure 11.— Dependence of the second scalar derivative-vorticity correlation, PyQ U)2 (43), on Taylor-microscale Reynolds number. 



(a) Dependence of the mixed scalar-velocity derivative correlation, F u g (40), on Taylor-microscale Reynolds number. 


Figure 1 2.- Scalar-velocity derivative correlations. 



(b) Dependence of the first scalar dissipation-strain correlation, F^jq 6 , (44), on Taylor-microscale Reynolds number. 


Figure 12.— Continued. 
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(c) Dependence of the second scalar dissipation-strain correlation, FyQ e2 (45), on Taylor-microscale Reynolds number. 
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Figure 12.— Concluded. 



(a) About 500 lines are plotted. A strong vortex structure extends most of the way across the box. The circle 
indicates where the vorticity of the structure is just below the threshold. From this perspective the vor- 
tex structure appears to be a tube. The arrow indicates where a two-dimensional slice is taken (figs. 15, 
16, and 18). 

Figure 13.— Three-dimensional perspective plots of the vorticity field from run F19 ( R ^ = 48.5). Short lines 
two mesh spacings long in the direction of the vorticity are plotted at the points of a 64 3 mesh if the 
vorticity is above a threshold. 
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(b) Same as (a), but rotated 90° about the vertical axis. In the middle the vortex structure continues to look 

like a tube. In the upper right it appears more twisted. 


Figure 13.— Continued. 
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(c) Same perspective as (a), but 2000 lines are plotted to show that the vortex structure does extend across 
the box. The circle indicates where two vortex structures might be merging into one. 

Figure 13.— Concluded. 
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(a) The same perspective as in figure 13(a). The scalar-gradient structure in the top center is aligned along 
with the vortex structure in figure 13 with the scalar gradient perpendicular to the vorticity. This scalar- 
gradient structure is a sheet seen edge on. Two scalar structures in the upper-left corner are sheets seen 
from the top down. 

Figure 14.— Three-dimensional perspective plots of the scalar-gradient field for the same time and run (FI 9) 
as figure 13. Short lines two mesh spacings long in the direction of the scalar gradient are plotted at the 
points of a 64 3 mesh if the scalar gradient is above a threshold. 500 lines are plotted. 
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(b) The same field as (a), but rotated 90° about the vertical axis. It has the same perspective as figure 13(b). 
The scalar-gradient structure in (a) which is aligned with the vorticity in figure 13(a) is the fan-like struc- 
ture in the upper center. The structures to the right and above the fan-like structure are the sheets in the 
upper-left corner in (a), but now seen edge on. 

Figure 14.— Concluded. 
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Figure 15 



Two-dimensional slice of vorticity as indicated by the arrow in figure 13(a). The structure in the 
center might be associated with vortex merging. 
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Figure 16 - Two-dimensional slice of the scalar gradient (figure 14) indicated by the arrow in figure 13(a). 
The structure in the center is partially wound around the vortex structure in figure 15 and shows the 
scalar gradient aligned perpendicular to the vorticity. 
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(a) The same perspective as in figures 13(a) and 14(a). The rate-of-strain structure in the top center is aligned 
along the vortex structure in figure 13 and the scalar-gradient structure in figure 14 with the compressive 
component of the rate of strain perpendicular to the vorticity and parallel to the scalar gradient. 

Figure 17.— Three-dimensional perspective plots of the compressive component of the rate-of-strain field for 
the same time and run (FI 9) as figures 13 and 14. Short lines two mesh-spacings long in the direction of 
the compressive component are plotted at the points of a 64 3 mesh if the total rate of strain is above a 
threshold. 1000 lines are plotted. 
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(b) The same field as (a), but rotated 90° about the vertical axis. It has the same perspective as figures 13(b) 
and 14(b). The rate-of-strain structure in (a) which is aligned with the vorticity and scalar gradient is the 
fan-like structure in the upper center. 


Figure 1 7.— Concluded. 
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Figure 18.- Two-dimensional slice of the compressive component of the rate of strain (fig. 17) indicated by 
the arrow in figure 13(a). The structure in the center shows the compressive component aligned perpen- 
dicular to the vorticity in figure 15 and parallel to the scalar gradient in figure 16. 
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